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An  Analytical  and  experimental  study  was  conducted  to  develop  a  computational  model 
adequate  for  stress  analyses  of  AS4/3501-6  composite  plates  containing  preformed 
1/4"  diameter  holes,  which  were  found  in  previous  and  parallel  experimental  studies 
to  possess  superior  open-hole  tension  and  compresssion  strengths  relative  to 
composites  with  drilled  holes  of  equal  diameter.  The  preformed-hole  composites 
contain  embedded  material  discontinuities  resulting  from  the  interpenetration  of 
adjacent  composite  layers  in  the  actual,  laminated  systems.  The  effort  reported 
addresses  the  initial  development  of  the  model  for  single  layers  containing  the 
material  discontinuities  present  in  laminates.  The  analytical  method  is  suitable 
for  adaptation  to  problems  involving  smart  meterials,  electronic  packaging,  micro¬ 
electromechanical  systems  and  joining/fastening. 

The  proposed  computational  model  is  based  on  spline  approximations  of  the  displace¬ 
ment  components.  Solutions  were  obtained  relative  to  curvilinear  coordinates  that 
mapped  the  preformed  hole  region  onto  a  rectangular  region  and  the  embedded  material 
boundaries  into  coordinate  lines.  Functions  capable  of  representing  point  discon¬ 
tinuities  in  the  strain  field  were  used  to  satisfy  the  equilibrium  conditions  at  the 
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interfaces  associated  with  the  embedded  material.  The  results  agree 
qualitatively  with  the  substantial,  experimentally  observed  increases  of 
the  tension  and  compression  strengths  of  plates  with  preformed  holes 
compared  to  those  with  drilled  holes. 


ABSTRACT 


Stress  analysis  in  the  vicinity  of  the  preformed  hole  with  the  fibers  smoothly 
modeled  around  the  hole  was  performed.  At  this  stage  the  stress  analysis  in  a  single  ply 
was  conducted.  Spline  approximation  of  displacements  was  employed.  The  solution  was 
obtained  in  curvilinear  coordinates  which  mapped  the  preformed  hole  region  onto  a  simple 
rectangular  region.  Two  predominant  fiber  orientations  were  considered:  0°  and  90° 
AS4/3501-6  plates,  both  containing  1/4"  diameter  preformed  hole  and  subjected  to 
unidirectional  tension  in  0°  direction.  For  the  0°  ply,  calculations  were  made  assuming  that 
the  matrix-rich  triangular-shape  gaps  left  by  the  fibers  surrounding  the  hole  are  filled  with 
90°  or  45°  fibers.  And  in  the  case  of  90°  plate,  it  was  assumed  that  the  gaps  are  filled  by  0° 
or  45°  fibers.  The  stresses  in  fiber  direction  and  in  the  direction  perpendicular  to  fibers  were 
compared  to  those  at  the  edge  of  the  drilled  hole.  At  least  40%  stress  reduction  in  the  0°  ply 
case  was  observed  for  the  stress  in  fiber  direction  at  the  edge  of  the  preformed  hole 
compared  to  the  same  stress  at  the  drilled  hole  edge.  The  work  was  supported  by  AFOSR 
task  #2303/DW  ,  "Biotechnology  for  Aerospace  Material  Requirements."  The  project 
engineer  was  Dr.  Richard  Hall,  WL/MLBM,  Wright-Patterson  AFB,  OH  45433-6533.  The 
contract  monitor  was  Dr.  Douglas  S.  Dudis,  WL/MLBP,  WPAFB. 

INTRODUCTION 

The  application  of  biological  ideas  to  improve  the  performance  of  composite 
materials  has  attracted  researchers  for  a  long  time.  Several  papers  were  devoted  recently  to 
the  strength  analysis  of  the  composites  containing  preformed  holes  with  continuous  fibers. 
In  the  paper  [1]  a  stream  line  function  is  utilized  to  describe  the  fiber  configuration  around 
the  moulded  hole.  A  glass/polyester  woven  roving  composite  was  considered. 
Unidirectional  lamina  with  moulded  hole  was  analyzed.  Unfortunately  the  calculated 
stresses  were  not  shown  in  the  principal  material  direction  and  didn't  allow  any  conclusions 
to  be  made  upon  the  strength  of  the  composites.  A  simple  equivalent  radius  model  for 
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estimation  of  the  compression  strength  of  the  woven  composite  with  the  moulded  hole  was 
proposed  in  [2].  Point  and  average  stress  criteria  were  applied  for  strength  prediction.  To 
obtain  the  stress  concentration  factor  at  the  moulded  hole  edge,  "Equivalent  Radius  Model" 
was  applied.  For  the  moulded  hole  of  a  given  radius,  an  equivalent  radius  of  a  drilled  hole 
was  obtained  so  that  the  compression  strength  for  these  two  cases  was  equal.  The 
composites  with  moulded  holes  were  found  to  have  higher  compression  strength,  therefore 
the  diameter  of  the  equivalent  drilled  hole  is  smaller  than  the  diameter  of  the  actual  moulded 
hole.  After  the  equivalent  diameter  is  found,  the  panel  with  the  moulded  hole  is  treated  as  a 
panel  with  a  drilled  hole  but  with  a  smaller  hole  diameter.  An  interesting  problem  is 
approached  in  [3]:  Can  the  tensile  capacity  and/or  the  buckling  capacities  of  a  plate  with  a 
central  hole  be  improved  by  using  the  curvilinear  fiber  format?  What  mechanisms  are 
responsible  for  this  improvement?  A  finite  element  approach  was  used.  In  every  element  a 
different  direction  of  principle  material  axis  was  considered.  The  results  of  stress 
calculations  were  compared  with  standard  quasi-isotropic  design.  It  was  shown  that  the 
curvilinear  fiber  format  has  the  potential  of  improving  the  load  carrying  capacity  of 
composites. 

A  common  disadvantage  of  the  models  considered  above  is  that  they  are  not  based 
on  considering  the  real ,  multidirectional  type  of  reinforcement  which  occurs  around  the 
preformed  holes.  Only  the  analysis  based  on  a  thorough  experimental  investigation  of  the 
structure  of  the  composite  with  the  preformed  hole  may  allow  to  understand  the  mechanism 
of  the  increased  or  reduced  strength  of  these  composites  depending  upon  the  applied 
loading.  Only  in  the  experimental  paper  [4]  was  the  multidirectional  reinforcement 
discussed  which  was  provided  by  the  expansion  of  adjacent  plies  into  voids,  or  matrix-rich 
regions  formed  by  the  fibers  surrounding  the  hole.  The  additional  structural  investigation 
and  photographs  provided  by  WL/MLBM  made  possible  the  development  of  the  model 
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proposed  in  the  present  report  and  used  for  the  stress  analysis  in  single  ply  composites  with 
moulded  in  holes. 

PROBLEM  FORMULATION 

Different  plies  of  the  de-plied  16-ply  AS4/3501-6  [90/45/-45/0]2s  specimen  with  a 
1/4"  preformed  hole  are  shown  schematically  in  Figure  1 .  Every  ply  can  be  divided  into 
two  regions.  Region  #1  is  the  region  with  approximately  constant  fiber  volume  fracture 
and  fiber  orientation  except  in  the  vicinity  of  the  hole.  According  to  [4]  a  4%  increase  of 
fiber  volume  fracture  in  the  local  preformed  hole  region  was  found.  The  continuous  fibers 
surrounding  the  hole  region  indicate  a  significant  local  change  of  fiber  orientation. 
However,  according  to  the  photographs  provided  by  WL/MLBM,  at  distances  more  than 
0.1  -  0.15'  inside  region  #1  from  the  gap  boundary  the  fiber  orientation  is  constant  for  the 
entire  ply.  The  region  #2  consists  of  two  triangular-shape  regions  which  do  not  contain  the 
fibers  from  region  #1.  A  sectioning  study  of  the  laminated  specimen,  performed  by 
WL/MLBM,  showed  that  these  regions  do  not  remain  pure  matrix  regions  in  the  laminate. 
They  are  filled  up  with  the  fibers  of  adjacent  plies  in  the  process  of  cure.  It  results  in  a 
multidirectional  in-plane  type  of  reinforcement  around  the  hole  which  is  responsible  for  the 
increase  of  the  preformed  hole  tensile  strength  as  compared  to  the  drilled  hole. 

In  the  present  report  the  stress  analysis  in  single  plies  with  preformed  holes  shown 
in  Figure  1  was  accomplished.  To  perform  a  satisfactory  laminate  analysis,  a  significant 
development  of  the  computational  procedure  was  needed  and  went  beyond  the  time 
limitations  defined  by  this  contract. 

METHOD  OF  SOLUTION 

Displacement  spline  approximation  in  combination  with  curvilinear  coordinate 
transformation  was  utilized  to  solve  the  problem.  Consider  an  orthotropic  rectangular  plate 
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of  length  -  L,  width  -  A  and  thickness  -  H  containing  a  circular  open  hole  (Figure  2a).  The 
gap  boundary  is  formed  by  straight  lines  emanating  from  line  y  =  yc  at  a  given  distance  D 
from  the  center  of  the  hole  and  ending  at  the  circular  hole  edge  tangential  to  it. 
Displacement  boundary  conditions  are  applied  at  the  edges  x  =  0,  L  so  that 
-  ux  (0,y,z)  =  ux  (L,y,z)  =  A. 

A  >  0  is  the  given  displacement  of  the  edges  generating  a  tensile  uniaxial  stress  oq 


Wf 

H  Aj0  l 


^XX  (0,y,z)  dydz 


The  following  curvilinear  transformation  of  the  x  and  y  coordinates  is  utilized: 

5 

x  =  Xc  +  ^  Fj(r)  X;(<j>), 
i=l 

5 

y  =  yc  +  XFi(r)Yi«to  (2) 

i=l 

-1  £r£  ls  0  <  «J>  <  2n 

The  curvilinear  coordinates  r  and  <))  are  generalized  coordinates  which  transform  the 
plate  with  the  open  hole  into  a  rectangular  region  [-  1,  1]  x  [0,  2tc].  It  is  very  important  that 
the  gap  boundary  in  these  new  coordinates  is  simply  a  coordinate  line  r=0.  The  open  hole 
boundary  is  the  coordinate  line  r=-l,  and  finally  the  coordinate  line  r=l  is  the  rectangular 
boundary  of  the  plate. 

The  functions  Fi  (r)  are  defined  by  the  following  simple  formulas: 
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Fi(r)  = 


fr2,~l<  r  <  o| 
(0,1  >r>()| 


F3(r)  = 


J-2r(r+l),-l<r<o\ 
2  F  \  0,  l>r>0  ) 

((r+1)2, -1<  r  <  0\ 

|  (r-1)2  1  >  r  >  0  / 

I  0,  -1<  r  <  0  \ 

F4(r)  = 

\2r(l-r),’  >r>Oj 

fo,  -1<  r  <  0\ 
Fs(r)  =  |r2,  1  >  r  >  oj 


(3) 


Graphics  of  these  functions  are  shown  on  Figure  2b.  The  functions  Xj(<())  and  Yi(<p)  are  as 
follows: 

Xi  (<J>)  =  d  cos  (<>  -  4>o) 

Y  i  (<j>)  =  d  sin  (0  -  <|)o) 


The  functions  X3  (<J>)  and  Y3  (<J>)  as  well  as  X5  (0),  Y5  (0)  are  expressed  in  terms 

of  cubic  splines,  providing  that  they  will  describe  the  gap  boundary  and  the  rectangular 

boundary  of  the  plate.  These  functions  are  shown  in  Figures  3a  and  3b.  The  functions 

X  2  5  Y2  and  X4 ,  Y4  are  not  independent  and  are  defined  as 

X2(0)  =  ^  (Xi  X3) 

Y2(0)  =  kY1+Y3) 

(4) 

X4(0  >  ~  2  x3  -  2  x, 
y4(0)=|y3-^-y1. 

One  can  obtain  that  the  relations  (4)  provide  that  the  derivatives 

9x  dy 

dr  dr 
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8F3 

are  continuous  for  any  r  and  0  values  ,  even  though  &  is  discontinuous  at  r  =  0.  The 
coordinate  lines  of  this  transformation  are  shown  in  Figures  4a  and  4b  for  the  0°  ply  and 
90°  ply  with  moulded  holes. 


STIFFNESS  PROPERTIES 

Let  us  consider  the  unidirectional  ply  material  constants  given  in  the  matrix  form 
C(123).  The  ply  properties  in  a  rotated  coordinate  system  xyz  can  be  obtained  according  to 
the  formula: 


C(xyz)  -  Te*(cc)  C(i23)  Te(a) 

Tg(a)  is  the  strain  vector  rotation  matrix  and  a  is  the  angle  between  the  x  -  axis  and  fiber 
direction.  This  matrix  is  given  in  reference  [5].  The  star  denotes  the  transpose  of  the 
matrix 


In  order  to  describe  the  stiffness  change  around  the  hole  region  due  to  the  local 
change  of  fiber  orientation,  the  following  assumptions  were  made.  The  Iocs’  increase  of 
the  fiber  volume  fraction  by  4%  in  the  vicinity  of  the  hole  was  neglected.  The  change  of 
the  fiber  orientation  was  described  as 


<J(r,<t>)  = 


9  -  ply  orientation,  r  >  0.2  \ 

e^  +  P(<j>)^,°>r>0.2) 


(5) 


In  formula  (5),  a(r,  <J>)  is  the  angle  between  the  principal  material  axis  direction  and 
the  x  -  axis  at  the  point  (r,  <j>),  0  is  the  ply  orientation.  According  to  (5),  a  is  a  constant 
value  for  r  >  0.2.  For  r  <  0.2  the  angle  a  is  smoothly  changing  from  0  to  (3  (0)  at  the  gap 
edge.  (3  (0)  is  the  declination  of  the  gap  boundary  to  the  x  -  axis  for  a  given  value  of  0.  It 
can  be  obtained  as 
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/dX 3\ 


(3(<}>)  =  arctg 


d<t> 


dY3 


d0  / 


The  coordinate  line  r  =  0.2  is  shown  in  Figure  3.  The  contour  defined  by  r  =  0.2  is 
approximately  at  a  distance  0.1-0.15"  from  the  edge  of  the  gap.  Thus  for  r  >  0.2  the  fiber 
orientation  is  constant  and  equal  to  ply  orientation  0. 

Considering  the  region  inside  the  gap  -  1  <  r  <  0,  we  assume  that  it  is  filled  with  the 
same  unidirectional  AS4/3501  but  having  the  fibers  oriented  in  the  direction  of  the  adjacent 
plies  in  the  laminate,  i.e.,  we  will  consider  the  following  cases: 
for  0  =  0' ,  a  =  90*  or  45°  if  - 1  <  r  <  0 

for  0  =  90°,  a  =  45°  or  0°  if  -l<r<0 


DISPLACEMENT  SPLINE  APPROXIMATION 

In  order  to  solve  the  complicated  problem  formulated  above,  we  are  using  the 
displacement  spline  approximation  approach  [6].  It  is  important  to  emphasize  that  by 
introducing  the  curvilinear  transformation  (2)  we  achieved  two  advantages:  first,  we 
transformed  the  complex  boundaries  of  the  plate  with  the  gap  into  simple  coordinate  lines 
r=  -  1,0, 1;  and  secondly  we  described  the  change  of  fiber  orientation  near  the  gap  using 
the  simple  relation  (5).  However  we  anticipate  great  difficulties  by  trying  to  describe  the 
orthotropic  stiffness  properties  in  the  curvilinear  coordinates.  Moreover  we  are  not 
interested  even  in  introducing  the  displacements  in  the  curvilinear  coordinates  because  these 
coordinates  are  just  general  curvilinear  coordinates  without  any  special  importance  from  the 
mechanical  point  of  view.  Therefore  the  unknowns  of  this  problem  will  still  be  the 
displacements  in  oxyz  coordinates:  ux  (x,  y,  z),  uy  (x,  y,  z)  and  uz  (x,  y,  z).  We  are 
going  to  write  the  strain  energy  expression  in  xyz  coordinates.  Then  we  will  replace  the 
integration  variables  x  and  y  by  r  and  0.  Let  us  introduce  the  following  vectors: 
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G(xyz)  —  (O’xxi  ®yy»  ®zz>  ®yz>  ^xz>  O'yx) 
£(xyz)  =  (£xx>  £yy«  £zz>  £yz>  £xz»  £yx) 


(6) 


The  standard  linear  relationship  between  displacements  and  small  strains  is  utilized. 

Hook's  law  can  be  written  in  the  oxyz  coordinates  as 

— »  -*• 

G(xyz)  =  C(xyz)  £(xyz)  •  (7) 

The  expression  for  the  strain  energy  is  as  follows: 


E  = 


£  C(xyz)  £  dxdy. 


where  A/Cl  denotes  the  rectangular  area  with  the  circular  hole.  The  stiffness  matrix  C(Xyz) 
is  changing  from  point  to  point  near  the  gap  in  accordance  with  (5).  The  strain  components 
in  this  formula  must  be  expressed  through  the  displacements.  Now  we  replace  the 
integration  variables  x  and  y  by  r  and  <{>  so  that 


E  = 


£*  Tg(a(r,<)>»  C(i23)  T£(a(r,<j>))  •  edet(J)  d<(> 


(8) 


To  complete  this  formulation  we  need  to  express  e  components  through  ux,  uy, 
Uz  in  r,  <J>,  z  coordinates.  The  Jacobian  can  be  represented  as: 

8x  dy 
dr  dr 
dx  dy 
,9$  dip/ 


In  order  to  express  the  e  components  through  ux  (r,  <j),  z),  uy  (r,  <|>,  z)  and 
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Uz  (r,  <)>,  z),  the  following  relations  can  be  used: 


(dux] 

(duA 

idlly  ^ 

0Uy 

(duz\ 

dx 

'  dr 

'  dx 

=  J» 

dr 

dx 

=jJ 

dr 

=  j1 

» 

> 

dux 

,3uy 

3uy| 

duZ| 

du2 

WJ 

Uy  j 

[tyf 

Uy/ 

i 

(9) 


Using  these  relations  we  can  write  (8)  in  terms  of  displacements  ux,  Uy,  uz  in  r,  <j>, 
z  coordinates.  Spline  approximation  of  displacements  upon  r,  <j)  and  z  is  used.  The  spline 
approximation  nodes  arc  introduced  in  r,  <}>  and  z  directions  so  that 

-  1  =  ro  <ri  <  . . .  <  rK  =  0  <  rK+  i  <• .  .<  rM  =  1, 

0  =  <f>0  <  4>1  <  .  •  •  <  $p  =  2jc 
0  =  zo<zi<...<zn=  H. 

Cubic  spline  functions  shown  in  Figure  5a,  5b  and  5c  arc  introduced  in  r,  <j>  and  z 

directions  respectively.  Note  that  the  set  of  functions  in  r  direction  incorporates  the 

discontinuity  of  slope  at  r  =  0,  which  is  the  boundary  between  two  materials.  In  the  4> 

direction  the  periodicity  conditions  are  required  at  <j>  =  0  and  <(>  =  27t.  The  displacement 

spline  approximation  can  be  written  as 

N+2  M+4  P 

Ux  (r,<M)  =  2  2  2  uijk  Zi(z)  Rj(r)  ®k(0) 

i=0  j=0  k=l 

N+2  M+4  P 

uy  (r,<}>,z)  =  X  2  2  vijk  Zi(z)  Rj(r)  Ok(<»  (10) 

i=0  j=0  k=l 

N+2  M+4  P 

UZ  (r,<J>,z)  =  222  WUk  Zi(z)  Rj(r)  ®k($) 
i=0  j=0  k=l 

The  equations  for  calculating  the  unknown  displacement  approximation  coefficients 
are  derived  by  substituting  (10)  into  (8)  and  requiring  the  first  variation  of  E  to  be  zero: 

8  (E)  =  0. 
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After  displacement  spline  approximation  coefficients  are  calculated,  the 
displacement  values  in  any  point  can  be  obtained  from  (10).  Then  the  strains  can  be 
calculated  using  (9)  and  the  stresses  can  be  obtained  using  Hook's  law  (7).  This  completes 
the  solution  procedure. 

NUMERICAL  RESULTS 

The  AS4/3501-6  unidirectional  ply  properties  utilized  in  calculations  are  given 

below: 

Ei  =  20  Msi  X  =  210  ksi 

E2  =  E3  =  1.7  Msi  Xd  =  210ksi 

G12  =  G13  =  0.9  Msi  Y  =  7.5  ksi 

G23  -  0.48  Msi  Yd  =  29  ksi 

V12  =  V13  =  0.32  S  =  13.5  ksi 

V23  =  0.55 

The  geometric  parameters  shown  in  Figure  2.a  were  obtained  by  measuring  the  specimens 
provided  by  WL/MLBM  and  found  to  be  equal  to 

d  =  1/4",  D  =  0.5",  L  =  5”,  A  =  2.5". 

Since  we  are  considering  a  single  ply  and  focusing  on  the  in-plane  stress  component 
calculation,  the  thickness  is  not  an  important  parameter.  H  =  0.04"  was  used  .  The 
displacement  applied  at  the  edges  x  =  0,L  was  equal  to  A  =  0.2  •  10'3.  It  resulted  in  the 

unidirectional  load  00  =  1595  psi  in  0°  ply.  Same  displacement  was  applied  to  the  90°  ply 
and  the  axial  stress  was  00  =  139  psi.  Since  ultimately  we  are  interested  in  stress  analysis 
in  the  laminate  we  will  normalize  all  stresses  in  both  plies  to  the  value  ao  =  1595  psi.  In 
Figure  6  the  stresses  at  the  gap  boundary  in  region  #1  are  shown  for  0°  ply  with  a  moulded 
hole.  Two  cases  are  considered:  when  the  gap  is  filled  by  the  unidirectional  composite 
with  fibers  in  90°  direction  and  in  45°  direction.  The  stresses  in  the  direction  of  fibers  at 
the  gap  boundary  are  compared  with  those  at  the  edge  of  the  drilled  hole  in  0°  ply.  The 
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angle  is  counted  counterclockwise  from  the  x-axis  direction.  Note  that  the  stresses  for  the 
moulded  hole  are  not  at  the  same  positions  as  in  the  drilled  hole  except  the  interval  from  78° 
to  112°  where  the  gap  edge  coincides  with  the  hole  edge.  However  even  there  the  fiber 
orientation  in  the  case  of  the  moulded  hole  is  continuously  changing  around  the  hole. 
Approximately  40%  stress  reduction  occurs  in  the  case  of  the  moulded  hole.  The 
importance  of  this  result  is  provided  by  the  fact  that  the  examined  stress  component 
controls  the  tensile  strength  of  the  0°  ply  with  the  hole.  The  reduction  of  the  stress 
concentration  indicates  that  one  reason  for  the  increased  tensile  strength  of  the  laminates 
with  moulded  holes  vs.  drilled  hole  is  more  efficient  stress  redistribution  within  the  0°  plies 
due  to  multidirectional  in-plane  reinforcement.  The  stresses  in  the  direction  perpendicular 
to  fibers  at  the  edge  of  the  hole  are  shown  in  Figure  7.  At  this  time  we  are  looking  at 
normal  stresses  in  the  direction  perpendicular  to  fibers  at  the  gap  boundary  ,  versus  same 
stress  component  at  the  edge  of  the  drilled  hole,  and  thereby  we  black  out  the  region  where 
the  gap  edge  coincides  with  the  hole  edge.  Higher  stresses  are  indicated  in  case  of  the 
moulded  hole  than  for  the  drilled  hole.  Practical  importance  of  this  conclusion  is  however 
questionable  because  it  may  be  a  mathematical  effect  due  to  singularity  at  the  point  where 
the  gap  edge  and  the  hole  edge  meet ,  forming  a  sharp  edge  in  region  #2.  The  results  for 
90°  ply  are  presented  in  Figures  8  and  9.  In  this  case  the  failure  is  obviously  controlled  by 
the  stresses  in  region  #1  in  the  direction  perpendicular  to  fiber  which  is  the  loading 
direction.  As  well  as  for  the  0°  ply  a  significant  reduction  of  the  critical  stress  component  at 
the  edge  of  the  moulded  hole  is  evident  as  compared  to  stress  at  the  edge  of  the  drilled  hole. 
Figure  8.  In  Figure  9  the  stress  in  fiber  direction  in  region  #2  is  shown  at  the  hole  edge. 
The  highest  stress  is  indicated  for  the  moulded  hole  with  the  gap  filled  with  0°  fibers. 
However  if  we  imagine  the  0°  and  90°  ply  with  in  a  laminate  and  compare  the  maximal 
Oi  i/ao  in  Figures  6  and  9  we  can  see  that  the  stress  in  the  90°  ply  with  moulded  hole  is 
much  smaller  than  in  0°  ply  and  is  not  a  factor  of  strength  reduction. 
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CONCLUSIONS 

1.  A  method  based  on  displacement  spline  approximation  and  curvilinear 
coordinate  transformation  was  proposed  and  developed  for  stress  analysis  in  single 
orthotropic  ply  with  moulded  hole.  The  region  of  the  gap  or  void  formed  by  fibers 
surrounding  the  hole  can  be  filled  with  the  unidirectional  composite  with  arbitrary  fiber 
orientation. 

2.  Stress  analysis  in  two  unidirectional  plies  with  fiber  orientation  of  0°  and  90° 
with  moulded  holes  under  uniaxial  tension  in  0°  direction  was  performed.  The  results  were 
compared  to  the  stresses  at  the  drilled  hole  edge  calculated  using  a  previously  developed 
program.  In  the  0°  ply  approximately  40%  of  the  reduction  of  the  failure  controlling  stress 
component  (stress  in  fiber  direction)  was  obtained  as  compared  to  the  drilled  hole.  For  90° 
ply  a  15%  reduction  of  the  failure  controlling  stress  component  (stress  in  direction 
perpendicular  to  fibers)  was  observed. 
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Figure  3.  Curvilinear  transformation  functions  in  circumferential 
direction 
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Figure  4.  The  curvilinear  transformation  in  the  0°-ply  (a)  and  90°-ply  (b).  The 
coordinate  lines  r  =  const. 


Spline  Functions  R^r)  tor 

7  Intervals  between  -1  and  0  and 
7  Intervals  between  0  and  1 


Spline  functions  <Dk($)  for 
17  intervals  between  0  and  2it 


Figure  5. 


Spline  functions  in  r  and  <}>  directions  :  (a)  and  (b)  respectively 
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Stresses  in  fiber  direction  in  0°  ply 


Figure  6.  Stress  in  Region  #1  for  0°  ply  at  the  gap  edge  with 
molded  hole  and  at  the  hole  edge  for  drilled  hole. 
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Stresses  in  the  direction  perpendicular  to 
fibers  in  0°  ply 


Figure  7.  Stress  in  the  Region  #2  of  0°  ply  at  the  hole 
edge  with  molded  hole  find  at  the  hole  edge  for 
drilled  hole. 
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Stresses  in  the  direction  perpendicular  to 

fibers  in  90°  ply 


—Gap  filled  with  0°  fibers 
-Drilled  hole 

—Gap  filled  with  45°  fibers 


Figure  8.  Stress  in  Region  #1  at  the  gap  edge  for  90 
ply  with  molded  hole  and  at  the  hole  edge 
for  drilled  hole. 
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Figure  9.  Stress  in  Region  #2  at  the  hole  edge  for  90° 
ply  with  molded  and  drilled  hole. 
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